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Introduction 

Financial modelling has been an area of extremely rapid growth in the past 30 years, and some 
extremely interesting mathematical challenges have emerged. One of the utmost importance for 
real-life applications to derivatives trading is that of calibration. Similar to common situations in 
many areas of physics and engineering, once a model has been suggested, its parameters have to 
be estimated using external data. In the case of derivative modelling, those data are the liquid 
(tradable) options, generally known as the "vanilla" products. It is well known since the pioneering 
work of Litzenberger and Breeden HI and its celebrated extension by Bruno Dupire [3] that the 
knowledge of market data such as the prices of vanilla options across all strikes and maturities is 
equivalent to the knowledge of the risk-neutral marginals of the underlying stock distribution, and 
moreover, that there is a unique one-dimensional driftless diffusion which recovers exactly such 
marginals. However, it has also been well-known for almost as many years that the evolution in time 
of the so-called "local volatility" is not stable, thereby leading researchers and financial engineers 
to look for a more robust, stochastic volatility type of modelling. In this paper, we consider the 
calibration problem for a generic stochastic volatility model: more precisely, we address the issue 
of calibrating to market data a generic model with a stochastic component and a local component for 
the volatility process. Such models are very useful in practice, since they offer both the flexibility 
and realistic dynamics of stochastic volatility models, and the exact calibration properties of local 
volatility models. In mathematical terms, the problem we consider is a non linear partial integro- 
differential equation for which we are able to prove short-time existence of classical solutions under 
suitable assumptions. The paper is organized as follows: Section 1 is devoted to the mathematical 
formulation of the problem. Section 2, to notations and statement of the main result. In Section 3, 
we recall some important technical results stemming from the general theory of parabolic PDE's. 
Section 4 contains the proof of the main result. Finally, Section 5 is a short conclusion. 



1 The Local and Stochastic Volatility model and its calibration 

The LSV model is an extension of the Dupire local volatility model. In the simplest situation - the 
two-dimensional case - the dynamics of the model are given by the following system of SDE's 

dS 

—± = a(t,S t )b(Y t )dBj +n t dt 
>j? 

dY t = a(tJ t )dBj + ^ t dt 



l 



1 f d 

-(-a 2 S 2 ( J b 2 pdy)) + -(rSq)+rq = 



Here, (S t ,t > 0) is the stock price process and (Y t ,t > 0) the stochastic component of the volatility. 
The function b simply transforms that factor into a proper volatility, a is the local volatility part of 
the model, choosing its value properly will enable us to calibrate the vanillas of the model, a is the 
volatility of the volatility factor and /j and £, are drift terms that may depend on the state variables 
and on time. B l and B 2 are standard brownian motions with correlation p. 

In order to fit the vanillas of this model, we write the Kolmogorov forward equation on the joint 
density p{t,S,y) of the couple {S t ,Y t ) 

I - w^ 2 " 2 ^ - m {pabaSp) - P l 2 a2p)) + i (M + rs {rSp) + r " = 

p(S,y,0) = 8(5 = So,y = y ) 
with (So,yo) the initial conditions. Taking q = f pdy the marginal of S, we get the equation 

dt dS 2 

Using Dupire's results from [3], we know that q has to solve the following equation in order to fit 
perfectly the vanillas of the market 

dq d 2 , 1 9 * . 3 , „ . 

where Op is Dupire's local volatility and contains the information about the vanillas we want 
to reproduce. We identify the terms in this last formula. This gives us the value of a 2 (t,S) = 

<5 2 D {t, S) j b i pdy =o|(f, S) /jxpdy - Eventually, the joint density that calibrates the smile of our model 
is solution of the nonlinear partial integro-differential equation 

+ ^($p) + ^(rSp)+rp = 

The rest of this paper is devoted to the study of a more general n-dimensional version of this 
equation. 



2 Generalized equation and notations 

Throughout this article, we denote by < t < T the time- variable and by x = (xi,X2, ■■■,x n ) E 
Q. C R" the n-dimensional space variable where Q. is an open subset with a sufficiently smooth 
boundary (we will precise this notion later). When we consider the equation from a financial 
point of view, the first-variable x\ stands for the spot and the last n — 1 for the volatility. Hence, 
we write S = x\, y = (x2, ■■■,x n ). We also denote by D T =]0, T[ x £1 the domain of definition 
and by B = {0} x a, B T = {T} x a and C T =]0, T[ x d£l the different parts of the boundary. 
Given the particular part played by the spot, we consider Q.s = {S E R/3y E R' !_1 , (S,y) E £1} and 
ws e ^ 5 ,^y = {ye w- l /(s,y) E n}. We are interested in the following equation: 

~ E ^-(P'7 a i«yP) + E^:(M+YP = on D T UB T (1) 

(J=2 ' J !=1 ' 



where (p ;j ) i<ij<n is a correlation matrix ie is symetric positive definite and verifies: pa = \ for all 
i and — | < p^- < I for i ^ j. We add the boundary condition p = *F on B U C T with *F constant on 
C T . We also let po denote the function po(t,S,y) = ^(S.y). The complexity of this equation stems 
from the following integral term: 

f a sp(t,S,X2,...,x n )dx2...dx n j a sp{t,S,y)dy 
fash 2 ( x 2 , ■ ■ ■ , x n )p{t, S, x 2 , . . . , x n ) dx 2 . . .dx n J^s b 2 (y)p(t, 5, y) dy 

(OQ) belongs to the class of nonlinear, parabolic and nonlocal equations. An interesting reference 
concerning that kind of equations is [7]. However, our case doesn't fall under the scope of that 
paper: the operator /(.) is not defined on Ct,(D T ). 
Let us now make a few remarks about our particular equation: 

1. in the case of an equation with no I(p) term, it becomes a classic linear equation of parabolic 
type. That kind of equation has been properly solved for quite some time now, see Il4l or ||5l. 

2. when b is constant, the problem is reduced to the previous remark. This observation is the key 
to our resolution method. First, we suppose that b does not vary too much and approximate 
the nonlocal term I(p) with a suitable constant. We then isolate the error made during this 
process in the second term and use a fixed point method to solve the new equation. 

3. in order to use this method and the results from flU, one has to assume that the coefficients of 
the equation belong to Holder spaces H k ' h > h ' 2 (we shall define them in the preliminaries). 

4. the question whether I is properly defined is natural. To answer it, we have to prove that 
f&sb (y)p(t,S,y)dy is bounded away from 0. To do it, we assume that b is non-negative and 
that *P the initial condition is strictly positive. By restricting ourselves to short times, we are 
sure that p is not too far from its initial condition and thus is strictly positive. 

5. from a financial viewpoint, it is natural to consider a domain Q. cylindrical with respect to 
the spot. However, since it may become very challenging to study a PDE on a domain with 
corners, we reduce our study to domains with a 5-section depending on S. 

The theorem we will prove requires the following assumptions on b and 

(Hi)fcec 1 0fr~ 1 ), 3(8i,52) Gl 2 ,0< 81 <b<8 2 onf"- 1 

(H2) V2 < i < n, |^| < b* on Q. y where b* is a constant we will choose later 

(H3) *F is strictly positive and in H 2,h,h / 2 . This gives us two results on po. First, po belongs to 
H 2 ' h ' h / 2 (D T ) and second 

< po = inf po < sup po = pb~ 

(H4) 0( V I / ) = on 35 in a sense described in the preliminaries 
Under the previous assumptions, we have the following result: 

Theorem 1. If the a,- belong to H 2,h,h / 2 (D T ), are positive and bounded away from by a stricly 
positive constant e, if the p\- are in H 1,h,h ^ 2 (D T ) and if J belongs to H 0,h,h ^ 2 (D T ), then, for b* small 
enough, there exists < T* < T and a solution of the equation ([/]) on D T UB T . 

The rest of the paper will be devoted to the proof of Theorem 1 above. 



3 Preliminaries 



In this section, we write d Xj and ^- without distinction. Let us start as in [4] and J51 with the 

following notion of distance d(P, Q) = [\x — x'\ 2 + \t — t'W 1 / 2 where P = (t,x) and Q = (t',x f ) belong 
to D T and \x\ is the norm of the n-dimensionnal vector x. Given such a metric d, we can define the 
concept of Holder continuity. For a function u, we write: 



sup\u 



ttD t , \ \u(P) — u(Q)\ , . n T . . n T rr r> T , \ 

H h («)= ™P Jn^/ \»\h =H0 + H h (") 



d t 1 p,q<ed t d{P,Q) 

H® 1 \u) < oo if and only if u is uniformly holder (exponent h) in D T . We denote by H 0,h > h / 2 (D T ) 

the set of all functions u for which \u\% T < °°. Now, if all the derivatives used in the equation exist, 
we write for k < 2: 

M& = Mf +z\d x u\f +-+n4u\f + \d t u\f 0) 

where the sums are taken over all the partial derivatives of the indicated order. We denote by 
H k ' h ' h ' 2 (D T ) the set of all functions u for which \u\® +h < °°. It is a Banach space and an algebra 
with the norm given by definition [3l Indeed, for all u,v in H k,h,h / 2 (D T ), we have: 

< HklhHth (4) 

We can now make the assumptions about D T more precise: for every point Q of C T , there exists an 
(n+l)-dimensional neighborhood V such that V fl C T can be represented, for some i (1 < i < n), in 
the formx, = r(t,x\, ...,Xi-\,Xi+i, ...,x n ) with r, d x r, d 2 x r, d t r Holder continuous (exponent h) and 
d 2 t r, df t r simply continous. 

We also have to consider functions \\r defined on BUC T . Such a function \\f is said to belong to 
H k,h,h/2 if there exists a ij/ in H k,h,h/2\pT^ such tnat vj/ = ^ on BUC T . We then define \\\f\ k+h = inf 

l^ljt+ft wnere me m f is taken with respect to all the ^'s in H k,h,h / 2 (D T ) which coincide with \|/ on 
BUC T . This process defines a norm on H k,h,h / 2 . 

The following results will be useful in the proof of our result. They concern the PDE: 

du n d 2 u n du 
Lu := - a ij( x ,t)-^- 5 \-J^bi(x,t)-^- + c(x,t)u = f(x,t) on D T UB T (5) 

hJ — * i— 1 

u = \\r on BUC T 

We shall need the assumptions: 

• the coefficients of the operator L belong to H Q,h,h l 2 (D T ), let be a bound on their norm 

n 

• for all (jc, t) in D T and for all ^eW\ £ ay (x, t)££j > K 2 \ £ | 2 (K 2 > 0) 

• V|/ G // 2 ^' A /2 and \f\D T < oo 

In addition, given the assumption about D T , if we consider a function \j/ e fj2,h,h/2^ f or an y ex t en . 
sion *P of \|/, 3^ is uniquely defined (by continuity) on the boundary dB of B, and the definition is 
independent of *P. We denote this function (on dB) by 3 f \|/. The other terms of L\)/ are also uniquely 
defined (by continuity) on dB. Thus, the quantity L\\f is well-defined on dB. 



Theorem 2. Under the previous assumptions and ifL\\f = f on dB, there exists a unique solution 
of the equation\5\ this solution belongs to H 2,h,h l 2 (D T ) and we have the Schauder inequality (with 
K H i depending only on K\, on K2, on h and on D T ) 

With < ^(M2 + fc + |/|?) (6) 
Furthermore, if\\f = 0, we can write a bound containing the time on the supremum of the solution 

Mo' < tKo\f\g (7) 

where Kq only depends on K\, on K2, on h and on £1. 

Proof: the first part of the result is classic, its proof can be found in [4J. As to the result with \\f = 
0, which is more original, one needs a result from [5] about volume potentials and representation of 
solutions of parabolic equations. It is the theorem (16.2) of section IV.16 we shall use. One reads 
that the solution of the equation [5] with \\f = can be written as 

u(x,t)= f dx f G(x,z,t,x)f(z,i)dz 
Jo Jo. 

where G is the Green's function for the operator L and verifies 

I |2 

\G{x,y,t,%)\ <K(t-x)-?exp(-K' lX ~ yl ) (8) 

t — x 

with K and K' two constants depending on the data of the problem. Using both these results, we 
get, for all t' < t and x E Q. 

\u(x,t')\<t\f\^ fdx! K{t-x)- n 2- l exp{-K l ^p^)dz<t\f\^K Q 
Jo Jo. t — x 

where Kq depends on K\, on K2, on h and on D T . 

4 Proof of Theorem d 

We are interested in the equation 

I - ^(P..a?/(p)rt - £^-( Pli a 1 a iv ^)p) 

- £ g-^-(Pyaiajp) + E^(PiP) +1P = on D T l>B T 

i,j=2 1 J !=1 ' 

and want to prove the 

Theorem. If the 0C ; belong to H 2,h,h / 2 (D T ), are positive and bounded away from by a stricly 
positive constant e, if the p\- are in 

H i,h,h/2( D T^ and if y b e i ongs to H°' h ' h / 2 (D T ), then, for b* small 
enough, there exists < T* < T and a solution of the equation (UJ) on D T UB T . 

Proof: 

The assumption (H2) gives us some control over the variations of b. Let us denote by b — b(yo) a 
strictly positive value taken by b (with yo £ Q.^ for some arbitrary S E Q5). We use the assumption 

on b to approximate the integral term I(p) with \jb}, the gap between those two quantities is 
quantified with the 



Lemma 1. There exists a constant Kb (depending only on h, n, 8\, 82, po and £1) and a polynomial 
function P strictly positive and increasing on such that \/p G H 2,h,h / 2 (D T ) verifying po < p, 
we have 

\I(P) - ^\2+h + \VKpj- \\th < b*K b P(\p\th). 

Remark. As a consequence of this lemma, we see thatVp G H 2 ' h > h / 2 (D T ) verifying p < p, I(p) 
belongs to H 2 > h h / 2 (D T ). 

We then write the equation as 

To solve this equation, we apply a fixed point method and use the lemma CD to get an upper bound 
on the second term. 

We take a real number x > Ipol^-u, an d t G and let X f x denote the set 

X< = {p G H 2 ^ 2 ^), \p\°' +h < x ^< p <pt+^, p = y onBUC T } 

The set X x clearly contains the function pq. We then consider the application M which takes a 
function u G X and sends it on v G H 2 > h,h l 2 (D T ) solution of the equation 

dv d 2 1 " d 2 1 " d 2 n d 

32 1 n d 2 1 

= ^(piia?(/(«) - ^)u) + ^^^(p li a l a i (^/I(^j--)u) (9) 

with the boundary condition v = ¥onBUC r . The existence of v is given by Theorem [21 Indeed, 
the coefficients of this equation belong to the appropriate spaces and because of (H4) the necessary 
condition 

oty = ^(pn«i(/w - + |^(p^^(Vw - lm 

on dB is verified. It remains to prove that this operator is elliptic: let (^;)i</<n be n real numbers 
and (t, S,y) G D T . We write f\ = ^ and f = a, for i > 2, we have 

£ p t jMt,S,y)fj(t,S,y)Z£j > KpfrffaS^g^Kpffg 

i,j=\ i=l (=1 

where the existence of K p is a consequence of p being a positive definite matrix. This proves the 
ellipticity of the operator, v exists and belongs to H 2,h,h / 2 (D T ). We now want to show that for 
suitable x and t, v belongs to X\ ie that 



For the first inequality, we apply [6] 



Mi'+h < ^(l¥l2+/ 1 + l^(Pii«i( / (")-^)") + E3^(Pi I -«i«Kv^-^'-'' 1 ^ 

T 1 n 1 / 

< ^2(|\K|f +/l + |piiaf(/( a )-^) a |^ + £|p h aia i ( v / ^- 7 -)M|?+J 



u 



h 



iL i=2 - 

< ^(M2+fc + (iPii«& 



h) 



< ^{\Mth + b*K'P{\u\t h )\u\th) (10) 

T n T 

where K' = (|pnOC 2 |f +/; + ^Ipi/CCia,!^)^ (we apply lemma Q] for the last line). 



7=2 

We remember that u belongs to X* x , thus \u\2 +h < x and then 

\v\th<K H 2(\y\? +h + b*K'P(x)x) 

Takings* = max{K H i(\^\f +h + 1), l/Jolf+z,) and ^* < ^P(l*)x* Ogives us |v|g^ < x*. 
It remains to prove that = < v < po + y . Let us write v = v — po- It is clear that v verifies 

32 j " 3 2 1 

O'v = O'p + ^p(p n al(I(u) - -£)u) + £^--(pi ! -aia ( -(^(«) - ^)") 

on Z) r U 5 r with v = on B U C T (here we use the fact that *F is constant on C ). 
We now apply the second part of Theorem [2l the inequality [7J to this function v 



3 2 1 " 3 2 1 

ft, o < ^l^o + ^(piiai(/(")-^)") + ^^(piiaiai(v^-^)")lo 



v 



£=2 

< f^odO'polf + |^(pna 2 (/( M ) - l)u) + f-l_( Pli aiai( - \)u)\% ) 



dS 2K 1V ' b 2 f^Sdxi 
< tK (\O'p \f + 1) (11) 

Taking T^od^'polo 7 + 1) = f ' we § et Ho' — T- Eventually, since v = po + v, the last inequality 

is proved and v belongs to Xj* . The application M maps Xj* into itself. 

Using this statement, we construct a bounded sequence (p n )neN °f functions belonging to Xj* 

• po has been previously defined 

• by induction, we write p n+ \ =M(p n ) 

By construction, we have Vn G N, \p„ < x* . Repeated applications of the Ascoli-Arzela theorem 
give us a function p E C 2 (D T ) limit in C 2 (D T ) of a subsequence of p n . Since 

I X t \-^U /) I 

r I OX; OX ; ^ ' ' OXiOXj^ ' ' 1 , , , / /, „T*-i 



We have 

I (x (x 1 I 

(u :y P + |,_V| r /2 .(»■'); (-•■')^')<.< 

And this computation being true for all the derivatives appearing in the norm H 2,h,h / 2 , we find that 

p £ £[2,h,h/2^ 

The last step of the proof is to take the limit in[9l The only result needed is I(p n ) — > Hp) . Since p n e 
X^* , the denominator is bounded away from 0. Two applications of the dominated convergence 
theorem give us the convergence we need. Thus, it is possible to take the limit in|9]which gives us 

— i—2 1 — i,J=2 ' J 1=1 ' 

= J^pnafc/cp) - ^)p) - l)p) 

That concludes the proof of the theorem, p is solution of our equation. 
Let us now prove lemma [IJ By definition, we have: 

\j( \-L\D T - i f pdy L\D T _ Jp(b 2 -b 2 )dy ]D T 

\W k 2^h-\ Jb 2 pdy ~ I k 2 Jb2pdy \2+h 



< 



Let us compute one after another the terms appearing in this norm (we remember that those func- 
tions only depend on t and S). Let (t,S) belong to ]0,T[ x Q5. We have 

I [ p(t,S,y)(b 2 -b 2 (y))dy\ < [ p(t,S,y)\b 2 -b 2 (y)\dy < \p\f [ \b 2 - b 2 (y)\dy 
JQS Jq$ JQ.S 

< IPlf [ , \b 2 (yo) - b 2 (y) \dy <b*{n-2)[ \y -y\dy\p\f 



< b*K\p 



T 



here and in the rest of the proof, K stands for some constant depending only on the data of the 
problem (Q., n, 81...) but not on p nor on b*. We get the last line from the following computation 
where y = (x 2 , .,x n ) and y = (V 2 , .,x' n ): 

n 

\b 2 (y )-b 2 (y) \ = \b 2 (x' 2 , .,xf n ) -b 2 (x 2 , .,x n )\ < £ \b 2 (x / 2 ,.,x' i ,Xi +1 ,.,x n ) -b 2 (x / 2 ,.,x' i _ 1 ,Xi,.,x n ) 

(=3 

< i\ d -p-\?\x' i - Xi \<(n-2)b*\y>-y\ 
(=3 0X1 



Now, let (t,S) and (t',S f ) belong to ]0, T[ x Q s . We compute 

| [ p(t^y)(b 2 (y )-b 2 (y))dy- [ P (t>,S>,y)(b 2 (y ) -b 2 (y))dy\ 

< f \p(t,S,y)-p(t',S\y)\\b 2 (y )-b 2 (y)\dy+ [ p{t^y)\b 2 {y Q ) -b 2 (y)\dy 
Jason* Jn§\n$ 

+ f^ s P(t,S,y)\b 2 (y )-b 2 (y)\dy 
<H° T (p)D((t,S),(t>,S>)) f \b 2 (y )-b 2 (y)\dy+\p\f f \b 2 (y ) - b 2 (y)\dy 

+ \P\o T f s> , \b 2 (yo)-b 2 (y)\dy 
<b*K\ P \f(D((t,S),(t',S'))+[ Jyo-y\dy+[. \yo-y\dy) 

By assumption on the boundary of our domain, it is possible to find a constant K depending only 
on D. such as V5,5' G Cls, Jqs^qs' \yo — y\dy < KD(S,S'). This gives us 



fp(b 2 -b 2 )dy\f <b*K\p\* 



And since | J p(b 2 -b 2 )dy\f +h = \ J p(b 2 - b 2 )dy\f + \ J^(b 2 - b 2 )dy\f + \ f%(b 2 - b 2 )dy\f + 
I / §s§(^ 2 — b 2 )dy\® T , we get from the previous computation 

\fp(b 2 -b 2 )dy\f +h <b*K\p\l h 

We now have to find a bound on I $tf-pdy 1 2+h~ Since P belongs to , we have I §d-pdy lo* — 
2 2 Now, let (t,S) and (t',S') belong to ]0,T[x Q S - We write 

1 1 \y l b 2 {y)p{S\t\y)d y -S sb 2 { y ) P {S,t,y)dy\ 

i i < y 



l f a sb\y)p(S,t,y)dy b 2 (y)p(S' ,f ,y)dy l ~ Snsb 2 {y)p{S,t,y)dy ^ b 2 (y)p(S',t',y)dy 

< KH°'( P )D((t,S),(t>,S>)) 

We used the same kind of arguments than earlier, K denotes here another constant depending on 8i, 
82, po and £1. This gives us I jj pc i y \h — ^0- + \p\% )• As for derivatives of jj pdy , f° r instance 
with respect to 5, we have 

\ d ( 1 \\D T -\ $ b2 % d y \D T < Wi 1 |„|D^2| n |D r 

l dS [ Wply~ )lh ~ l ~(fb 2 pdy) 2lh + > 

The same kind of computation is true for the derivative of second order 

|9 f !b 2 jdy T 2{fb 2 %dy) 2 fb 2 ^dy T ^ T T T T 

l dS { ~(fb 2 pdy) 2)lh "I (Jb 2 pdyf ~(Jb 2 pdy) 2lh + ^ } + + ^ ) I^W 

Eventually, we get 

I J^l£* < + \P\tk + {\p\tkf + (\p\f +h f + (\p\f +h ) 4 + (\p\f +h ) 5 ) 



Combining this result with the previous computations, we find 



\I(P) 




<b*KP(\p\? +h ) 




b 



, we find the same kind of results for the second term involved in the lemma. This con- 



cludes the proof. 
5 Conclusion 

In this paper, we have shown that the equation driving the calibration problem for local and stochas- 
tic volatility models is well-posed in the case of suitably regularized initial conditions. It is however 
clear that the solution of the full Kolmogorov equation with Dirac initial condition does not obtain 
as a consequence of Theorem 1 : possible extensions of our results towards this direction are cur- 
rently being explored. Let us also mention that a generalization of Theorem 1 to multidimensional 
correlation calibration have already been investigated and will be presented in (HI 
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